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Abstract 

This paper is a theoretical and a numerical investigation of the stability of a tip-cantilever system 
used in Non-Contact Atomic Force Microscopy (NC-AFM) when it oscillates close to a surface. 
No additional dissipative force is considered. The theoretical approach is based on a variationnal 
method exploiting a coarse grained operation that gives the temporal dependence of the nonlinear 
coupled equations of motion in amplitude and phase of the oscillator. Stability criterions for the 
resonance peak are deduced and predict a stable behavior of the oscillator in the vicinity of the 
resonance. The numerical approach is based on results obtained with a virtual NC-AFM developped 
in our group. The effect of the size of the stable domain in phase is investigated. These results 
are in particularly good agreement with the theoretical predictions. Also they show the influence 
of the phase shifter in the feedback loop and the way it can affect the damping signal. 

keywords : NC-AFM, Variational principle. Stability, Virtual machine. Phase shifter. Damping 
variations. 
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I. INTRODUCTION 



In recent years, the use of the Non-Contact Atomic Force Microscopy (NC-AFM) mode 
has shown that contrasts at the atomic scale could be achieved on semiconductors and insu- 
lators surface a^i^i^i^i^ . Experimental and theoretical features dedicated to the description of 
this dynamic mode have been widely discussed in previous papers^'^'^'SiiSiiiJ^. In particular, 
it was shown that the high sensitivity of the oscillating tip-cantilever system (OTCS) was 
based on the value of the quality factor and on its nonlinear dynamics in the vicinity of the 
surfacoiSiH. Current considerations of the authors focus on the origin of the increase of the 
damping signal when the tip comes close to the surface, in the range of a few angstroms. 
Some claim that the origin of this apparent increase could be due to the hysteretic behavior 
of the 0TCSi2»^. These interpretations implicitly rise the question of the stability of the 
OTCS when it is at the proximity of the surface. 

The aim of this paper is to show from a theoretical and a numerical point of view that the 
nonlinear dynamics of the OTCS leads to various stability domains of its resonance peak that 
may help to understand the reason why the NC-AFM mode, while being so sensitive, keeps, 
in most of stable behavior. In other words, this work is an attempt to show that, 

if no additional dissipative force is considered between the tip and the surface, an apparent 
increase of the damping signal cannot be the consequence of the nonlinear behavior of the 
OTCS. 

The paper is organized as follow. The first part is dedicated to a description of the non- 
linear behavior of the OTCS at the proximity of the surface. To do so, a specific theoretical 
frame based on a variational method using a coarse grained operation has been developed. 
This gives the explicit temporal dependance of the OTCS equations of motion. These equa- 
tions are the basis to analyze the stability of the stationary state. A large part of this work is 
detailed in ref.^^^. Experimentally, as the NC-AFM mode requires the use of a feedback loop 
to maintain constant the amplitude at the resonance frequency, which in turn requires to 
maintain constant the phase of the oscillator around its resonance value, e.g. —ti/2 rd, the 
phase variations of the OTCS will be extracted and discussed. The second part of the paper 
deals with numerical results obtained with the virtual NC-AFM^^ which is very similar to 
the experimental machine. These results show unambiguously the contribution of the phase 
shifter in the feedback loop and the way it can lead to damping variations. 
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II. THEORETICAL APPROACH OF THE NC-AFM 



The present section is divided into three parts. The first one details the specific theoretical 
frame for the obtention of the equations of motion in amplitude and phase of the OTCS. 
A coarse-grained method gives the equation describing the time evolution of the stationary 
state of the OTCS as a function of the coupling term between the tip and the surface. 
The second part is a description of the distortion of the resonance peak as a function of the 
distance. This part was detailed in ref.— so that only the main results are given. The results 
provide the basis of the discussion about the stability of the branches which is detailed in 
the third part. 

A. Theoretical frame 

We search a solution to the temporal evolution of the OTCS by using a variational 
solution based on the principle of least action. Even though this approach exploits the same 
physical concepts than the one which had led to the coupled equations in amplitude and 
phase of the stationary state of the OTCSi^ii^, it appears to be more general since here, the 
temporal dependance is explicitly obtained. We start from the definition of the action of 
the OTCS coupled to an interaction potential : 

S= [' C{z,z,t)dt, (1) 

Jta 

where C is the Lagrangian of the system and z (t) the position of the tip with time^^ : 



C{z,z,t) = T-V + W 



(2) 



-kcZ {tf - Z (t) Texc COS iut) + Vint [z if)] 



m luq 



z{t)z{t) 



ujq, Q, m* and kc = m*ujQ are respectively the resonance pulsation, quality factor, effective 
mass and cantilever' stiffness of the OTCS. J^exc and uj are the external drive force and drive 
pulsation. Due to the large quality factor, we assume that a typical temporal solution is on 
the form : 

z{t) = A it) cos [ut + if it)] , (3) 
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where A (t) and (p (t) are assumed to be slowly varying functions with time compared to the 
period T = lix juj. The underlined variables of z (t) in equ|21 e.g. A(t), ip (t), A (t) and (t), 
are calculated along the physical path, thus they are not varied into the calculations^^. 

To describe the interaction between the tip and the surface, the attractive coupling force 
is assumed to derive from a sphere-plane interaction involving the disperse part of the Van 
der Waals potentiaU^ : 

Vint [Z it)] = , , (4) 

^ntl K Ji 6[D-Z (t)] ^ ^ 

H, R and D are the Hamaker constant, the tip's apex radius and the distance between the 
surface and the equilibrium position of the OTCS. 

The equations of motion in amplitude and phase of the OTCS are obtained by considering 
the following coarse-grained operation. Let's assume a long duration At = th — ta with 
At ^ T and calculate the action as a sum of small pieces of duration T : 



.{n+l)T / . .(n+l)T \ 

S = J2 C{z,z,t)dt = J2(7f; C{z,z,t)dt]T = J2^eT (5) 

JnT „ y JnT J 



Ce is the mean Lagrangian during one period and appears as an effective Lagrangian for 
a large time scale compared to the period. Owing to the quasi-stationary behavior of the 
amplitude and the phase over the period, the effective Lagrangian is calculated by keeping 
them constant during the integration. The calculations give : 

£e(AA<^,0) = —[A^ + A^ {u + v3^)J - ^ + ^ T ^^nt[zit)]dt 

J 

AAcos [if — ip) — AA + (p) sin {ip — p) (6) 

Note that the effective Lagrangian is now a function of the new generalized coordinates A, p 
and their associated generalized velocities A, p>. At this point, remembering that the period 
is small regardless to At = th — ta during which the total action is evaluated, the continuous 
expression of the action is : 

S= r cjA,A,ip,ip)dT, (7) 



where the measure dr is such that T ^ dr <C At. 

Applying the principle of least action 6S = to the functional £e, we obtain the Euler- 
Lagrange equations for the effective Lagrangian. Thus, the amplitude and phase equations 



of motion of the OTCS coupled to the surface through an interaction involving the disperse 
part of the Van der Waals potential are obtained : 



r, .,2 ,1 a cos (09) an 



Q Q 3{d^-a^f^ /ON 
/2d IV sin(v;) ' ^""^ 



\ a Q J aQ 
In equIHl specific notations were used, d = D/Aq is the reduced distance between the 
location of the surface and the equilibrium position of the OTCS normalized to the resonance 
amplitude Aq = QJ^exc/K, o = ^Mo is the reduced amplitude, u = uj/uq is the reduced drive 
frequency normalized to the resonance frequency of the free OTCS and Ha = HR/ (^c^o) 
the dimensionless parameter that characterizes the strength of the interaction. 



B. Resonance frequency shift 

The equations of motion of the stationary solutions a and (p are obtained by setting 
d = = and d = = in equIHl and lead to two coupled equations of the sine and 
of the phase of the OTCS previously calculatec^^ : 



cosme 



cos ((/?) = Qa{l — 



aQKa 



3 (d^ — a 



(9) 



sin {(p) = —ua 

Solving equO gives the relationship between the frequency and the amplitude at a given 
distance dM : 



u± [a) 



It 



\ 



1 - 4g2 1 



(10) 



The signs plus and minus are deduced from the sign of cos (ip) and correspond to values of 
the phase ranging from to —90° cos (p) > 0) or from —90° to —180° cos (p) < 0), 
in agreement with the sign convention of the phase in equOI From equEl is calculated the 
resonance peak at any reduced distance for a given strength of the sphere-surface interaction 
and equini give the associated phase variations. The two branches define the distortion of 
the resonance peak as a function of d. U- gives the evolution of the resonance peak for 
frequency values below the resonance one and u+ for frequency values above the resonance. 
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In fig is given the distortion of tlie resonance peak vs. the reduced distance d. For large 
values of ci, e.g. when the surface is far from the OTCS, the nonlinear effects are negligible 
and the peak keeps a well-defined Lorentzian shape (see equ]l(JI with = 0). When the 
OTCS is approached towards the surface, because the interaction is attractive, the resonance 
peak starts to distort towards the low frequencies. The distortion of the peak increases as 
d decreases. On the figlH the branches that are supposed unstable are shown with dashed 
lines. 

Using equHni the resonance frequency shift as a function of the distance d is obtained by 
setting a = 1. This former condition ensures the required condition for the NC-AFM mode. 
Thus, the normalized frequency shift, {u — uq) jv^^ is given by m — 1— : 



The frequency shift given by equlTTI vs. d can be deduced from the figd (see the arrows on 
the figure). Following the previous discussion about the stability of the different parts of 
the resonance peak during the distortion, since the measure is performed as a function of d 
with a = 1, no bistable behavior can be observed. 

However, in the vicinity of the resonance the branches m+ and m_ become very close 
(see for instance figH with cis = 1.012). Therefore, even with an oscillation amplitude 
kept constant, question rises about the ability of the OTCS to remain on the same branch. 
Qualitatively, one may expect that around a = 1, the branch U- is unstable and is 
stable (see fig^). If this is true, any small fluctuation of the oscillation amplitude might 
produce a jump from one branch to the other one as discussed in refs .^° '^ ^ '. Since the branch 
u_ seems to be unstable, a jump to this branch should lead to an abrupt decrease of the 
amplitude, which in turn might produce an apparent abrupt increase of the damping signal 
as a consequence of the hysteretic behavior. Because such a jump is, in most of cases, never 
observed, it becomes useful to determine more accurately the stability of the two branches. 

C. Stability criterions 

The stability of the branches u± of the resonance peak is obtained from equations of mo- 
tion of the OTCS (see equIH}. By linearizing these equations around the stationary solution 




(11) 
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(now identified by the index "s" ) and using classical considerations of the linear theory, one 
gets the stability criterions of the branches u± . The stability criterions can be expressed 
from the derivatives das/du± of the branches and reduced to the simple expressioi>i^: 

' dd 

-j-^ > and cos (ips) > dsj (2(5) (0 

or (12) 
< and cos {^p^ < a^/ (2Q) (ii) 

K du 

The figsI2](^aj and (b) show the distortion of the resonance peak and of the associated phase 
curve, respectively. The figsOlf^aj and|ni(^&j are zooms on the region a of the figs|21(^aj and 
(b), respectively. 

For the branch da^/du^ being always negative and the associated value of the phase 
being always defined beyond — 7r/2 (see section lllB|) . the criterion {ii) implies that u+ is 
always stable, whatever the value of a^. 

For M_, the sign of the derivative changes twice. For this branch, the phase is always 
defined above — 7r/2. Therefore on the lower part of the branch (small a), dag/du^ > and 
the criterion (i) indicates that the branch is locally stable. When das/du- becomes negative 
(see figElfaj), because ips > — 7r/2 (see fig^b)), the criterion (i) is no more filled. As a 
consequence, m_ is locally unstable and the instability is precisely located where the infinite 
tangent appears. On the upper part of the resonance peak, the curvature of U- changes 
again and dag/du^ > (see figOir'^y')) implying that it is again a locally stable domain. 
Thus the branches U- and if- exhibit two stable domains and one unstable. 

Note also that the resonance condition is deduced from das/du = which implies 
cos (ips) = as/ (2(5), or equivalently u_ = u^ or again yj. = ip+. This equality is the 
usual resonance condition of a free harmonic oscillator. If = 1, e.g. without any coupling, 
the resonance phase is therefore ips = arccos[l/ i2Q)]. For the OTCS we used, Q ~ 500, 
and so (fs — — vr/2. But taking into account the fact that the coupling only slightly modifies 
the value of the resonance amplitude, ~ 1.0013 (see fig|2l(^aj), we still obtain (fs — —7^/2 
so that we can consider that the nonlinear resonance is always given by the relationship 
(fs = -vr/2. 

Therefore the theoretical approach foresees that u+ is always stable but that also a small 
domain of m_ around the resonance value remains stable. If the resonance value would have 
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been located at the point where dag/du- is infinite, an infinitely small fiuctuation would 
have been able to generate an abrupt increase of the damping signal as discussed previously 
and suggested in ref.~, or more recently in ref.^^. Experimentally, an electronic feedback 
loop keeps constant the amplitude of the OTCS so that its phase is located around — 7r/2 
(see section below). As a consequence, question rises about the size of the stable domain 
in phase around —7i/2. If any fiuctuation around —tt/2 makes the phase going beyond the 
stable domain, the OTCS behavior becomes unstable. For Q = 500, the size of the stable 
domain is of about 2.6.10^^ rd (see figEir^j) whereas it's reduced to 2.6.10^^ rd for Q = 5000 
(data not shown). Thus, if the electronic loop is able to control the phase locking with a 
better accuracy than 2.6.10"^ rd, the OTCS will be locked on a stable domain. 

Therefore, if the setpoint of the oscillator is properly located at the ~tt/2 value throughout 
an experiment, this value corresponds to a stable domain and consequently will neither give 
rise to amplitude or phase variations nor to damping variations. 

III. VIRTUAL NC-AFM RESULTS 

In a recent paper, we have described a virtual NC-AFM machine built using the Matlab 
language and the Simulink toolbox^^. This machine is very similar to our own experimental 
hybrid machine built with Digital Instruments^ and Omicron^^ blocks. The virtual ma- 
chine has been extensively used to study the frequency shift and the damping signal in the 
approach- retract mode. Two types of situations have been investigated : i— the first one 
corresponds to the case where no dissipative force is introduced in the tip-surface interaction, 
ii— the second one deals with dissipative forces. In both cases an attractive sphere-plane 
Van der Waals interaction is taken into account. 

In spite of previous results that have already shown that the damping signal could be 
considered as a constant when no dissipative force was introducedAS, here we want to inves- 
tigate with the virtual machine the stability of the OTCS by looking accurately at its phase 
variations within the electronic feedback loop that maintains constant the amplitude of the 
oscillations and compare these results with the theoretical predictions. To do so, we still do 
not consider any additional dissipative force. 

The theoretical results have led to the conclusion that, provided that the OTCS phase is 
in the close vicinity of — 7r/2, this setpoint corresponds to a stable branch. As a consequence, 
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the questions are : i— What is the part of the feedback loop that controls the size of the 
phase stable domain of the OTCS around -7r/2 rd ? and ii— Would it be possible to change 
the parameters of this element in order to change the size of the phase stable domain and 
thus observe phase variations? As a consequence of these changes, variations of the damping 
signal should also be observed. 



A. The phase shifter of the feedback loop 

In figEl is drawn a very simplified schematic diagram of the feedback loop of the NC- 
AFM (for more details, see ref.— ). Usually, the phase (tu) of the phase shifter transfer 
function is adjusted to — 37r/2 so that the loop oscillates at uq which is the free resonance 
frequency of the cantilever, corresponding to a tip-surface distance D ^ oo. We recall that 
the oscillations of the loop are ruled by the relation : 

(j){uj) + (p{uj) = 0±2mT, (13) 

where n is an integer and (uj) is the phase difference between the oscillations and the 
excitation of the cantilever. If the setpoint is fixed to the resonance frequency, then ip {uq) = 
—7i/2. The phase adjustment in the Omicron electronics is obtained by changing the bias 
of varicap diodes^S. The phase shifter transfer function in terms of the p Laplace variable 
can be written as H (p) = I , the time constant r being adjusted by the user such 



that, at the resonance : 



1 + Tp 



(f) {ujq) = —4 arctan {tuq) = — — (14) 

When the tip-surface distance D is reduced, due to the coupling, uq decreases. As a 
consequence, (f){uJo) and (f {ujq) and are no more equal to — 37r/2 and — 7r/2 respectively. 
According to equHSl the variation of (tu) is governed by the one of (t>{uj). Assuming a 
small variation around the resonance frequency Au = u — ujq, one gets : 

(uj) ^ -37r/2 - (15) 

As D decreases, Au is negative. Therefore (f){uj) becomes larger than — 37r/2 and ip{uj) 
smaller than — 7r/2. The decrease of <f{uj), f (uj) < — vr/2, means that the phase of the 



9 



OTCS follows the phase branch associated to u+, ip+, which is always stable (see fig|21f6j). 
Thus the loop is always stable. Moreover, the hypothesis implying that {uj) keeps a 
value close to — 7r/2 is a very good assumption. To proof that, let us consider for instance 
uq = 150 kHz, which is a reasonable value for a cantilever. Therefore r = 2.56.10"^ s (see 
equ ll4j) . Assuming now a large frequency shift, Au = —200 Hz, we get A0 (u) = 1.9.10^^ rd 
and therefore A(f{uj) = —1.9.10"^ rd. In spite of the rough assumption of a first order 
expansion of the phase of the phase shifter, the typical phase variations of the OTCS around 
the nonlinear resonance are less than 2.10^^ rd. This implies that the machine properly 
follows the nonlinear resonance, even when large frequency shifts are considered. 

The curves [a] in figsEl and El show the phase ip (u) and the damping signal vs. the 
distance D, respectively. As expected, the variations are very weak. 



B. "Controlled" damping variations 

If we want to observe the phase instability predicted by the theoretical calculations, 
the phase shifter transfer function should have been on the form d(f){uj) /duo > Q around 

ujq. a possible expression of such a transfer function would be : H {p) = 



1 + rp 



1 — rp^ 

Experimentally, this form is not feasible and even if it were, the loop would become unstable 
and therefore no stationary state could be reached. The reason is that the inverse Laplace 
transform of 1/ (1 — rp) varies as e*/'^ which diverges as t ^ oo. 

In the virtual machine, it is possible to implement a phase shifter with a slope d(j) {u) /duj 
larger than the phase shifter built by Omicron. We have retained the following transfer 
function which is easy to do with electronic components : 



H{p) 



VI 

9 , ^1 , 2 



(16) 



The parameters uJi and Qi may be adjusted to obtain for instance (I){ujo) = — 37r/2. The 
phase of the transfer function is then : 

0(.) = -2arctan(^-^^|^) (17) 
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For a small frequency shift, dcp (u) / du ^ —AQi/uq. Keeping the same values than previously 
uq = 150 kHz and Au = —200 Hz and assuming Qi = 50, we now obtain a change {uj) of 
about 0.26 rd. Consequently, the change Aip (u) becomes larger (see the curve [b] in figE)) 
and we now observe an increase of the damping signal as shown in figEl curve [b]. 

The previous examples are pedagogical cases for which an arbitrary large value of the 
slope of the phase of the phase shifter was considered. The ideal phase shifter should 
maintain the phase (p i^o) at — 37r/2 so that the frequency of the loop remains equal to the 
resonance frequency of the cantilever. This is not possible in practice, however it is clear 
that the solution retained by Omicron is very close to the ideal case {uiq) = — 37r/2 because 
d(j) (u) jdui is very weak. 

IV. CONCLUSION 

A variational method based on a coarse-grained operation has been used to investigate 
in details the stability of an oscillating tip-cantilever system near a surface. The tip-surface 
interaction is described by Van der Waals forces. Results show that the resonance peak of 
the oscillator can be described from two branches. The first one, named u+, corresponds to 
frequencies larger than the resonance. Stability criterions deduced foresee that it is always 
stable. The second one, m_, may be decomposed into three domains : two are stable and one 
is unstable. The second stable domain of u_ is small and is defined at the upper extremity 
of the resonance peak. The phase at the resonance ip {ujq) = —tt/2 is at the overlap of the 
u+ and of this former second stable domain of m_, thus the setpoint ip {uq) = — 7r/2 belongs 
to a stable zone. 

This result is of great importance to understand the stability in NC-AFM. In this tech- 
nique, the phase of the cantilever is adjusted to — 7r/2 within an electronic feedback loop 
as the tip-surface distance D is infinite. In the approach mode, the frequency of the loop 
decreases, consequently the phase becomes smaller than — 7r/2 because the slope d(f){uj) /duj 
of the phase of the phase shifter transfer function is always negative. Thus the oscillator 
always "slides" along u+ and the system is unconditionally stable. This is what is usually 
observed experimentally and confirmed by the results of the virtual NC-AFM we have built. 
Because the slope del) {uj) /doj and the frequency shift are very weak, we may consider that 
the phase ip {ujq) of the oscillator is always very close to — 7r/2, typical variations being less 
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than 2.10 ^ rd. Consequently, the damping signal keeps constant if no dissipative force is 
introduced in the tip-surface interaction. 
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FIG. 1: Distortion of the resonance peak computed from equ llOl for three values of the distance, 
di = 2, d2 = 1.11 and d^ = 1.012. The numerical parameters are = 20nm, Q = 400 and 
Ka = 8.10"'^. For an attractive coupling, the peak is more and more distorted towards the low 



frequencies as d is reduced, e.g. the surface is approached. For each value of d, the unstable 



domains of U- are shown with dashed lines. The arrows indicate the resonance frequency shift vs. 



d. 
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FIG. 2: (a)- Distortion of the resonance peak computed from eau llUI The numerical parameters 
are d = 1.05, Aq = Wnm, Q = 500 and Ka = 2.5.10"^. The stabihty criterions foresee that 
ii+ is always stable whereas U- exhibits two stable domains (continuous lines) and one unstable 
(dashed lines). The domains are separated by the spots where the derivative da/du- is infinite. 
(h)- Distortion of the phase curve computed from equsE] associated to the resonance peak. As a 
consequence of the stability of u+, 994. is always stable whereas 93_ exhibits two stable domains 
and one unstable. 
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FIG. 3: (a)- Zoom in the region a of the resonance peak. As da/du- becomes infinite again 
criterions define a new domain of U- which is stable. The resonance is located where u+ = 
(b)- Zoom in the region a of the phase curve. The resonance is located at —tt/2 where (/?_(_ = 
and belongs to a stable domain. 
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FIG. 4: Schematic diagram of the feedback loop used in the virtual NC-AFM which is very similar 
to the one of the experimental machine. 
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FIG. 5: Variations of the phase of the OTCS (p{uj) within the feedback loop vs. the distance 
D computed from the virtual NC-AFM. The numerical parameters are : resonance amplitude 
Aq = 15 nm, spring constant kc = 40 N.m^^, quality factor Q = 5000, tip's radius i? = 10 nm 
and Hamaker constant H = 2.10"^^ J. Curve [a] : The phase of the transfer function H (p) 

of the phase shifter is the one given by equll4| e.g. is similar to the experimental machine. As 
D decreases, ip (lo) becomes weakly smaller than — 7r/2 (less than 2.10^^ rd), therefore follows the 
stable branch The machine follows accurately the setpoint, which is always stable, even when 
the tip is in the very close vicinity of the surface. The associated damping variation nearly no 
varies (see figEI)- Curve [b] : (j){uj) is the phase of H (p) whose expression is given by eau ll61 
Around —tt/2, the slope d(j) (lo) /dw is larger than in case [a] so that ip [lo) decreases more quickly. 
As a consequence, the damping signal increases. 



17 



L44 - 
£ 1-^3 - 




1.39 I ' ' ' 1 

0.2 0.4 0.6 0.8 1 

Tip -surface distance (nm) 

FIG. 6: Variations of the damping signal vs. the distance D. Curve [a] : No damping variation is 
observed if the phase of the virtual NC-AFM phase shifter is similar to the one of the experimental 
machine. In the very close vicinity of the surface, as (p{ujq) ^ — vr/2, e.g. the amplitude of the 
oscillations slightly decreases, a weak increase is observed. Curve [b] : As (p (loq) varies more 
quickly than in curve [a] due to the different expression of H (p), a larger increase of the damping 
is obtained. 
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